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ABSTRACT 


The  dynamic  lift  and  drag  response  of  a  large  8  ft  (2.43  m)  span, 
aspect  ratio  6  hydrofoil  operating  in  regular  waves  was  experimentally 
determined  for  a  range  of  chord  Froude  numbers  varying  from  1.22  to  3.66. 

The  amplitude  of  the  fundamental  frequency  lift  response  was  found  to 
increase  significantly  with  decreasing  chord  Froude  number,  while  the 
second  harmonic  lift  response  was  zero.  The  amplitude  of  the  first  and 
second  harmonic  drag  response  also  varied  with  chord  Froude  number.  The 
first  harmonic  drag  response  was  a  maximum  at  the  lowest  chord  Froude 
number  of  1.22,  and  the  shallowest  depth  chord  ratio  of  0.5.  Several 
theoretical  predictions  of  the  lift  amplitude  response  are  compared  with 
the  experimental  results. 

ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  was  carried  out  for  the  Ship  Performance  and 
Hydromechanics  Block  Program,  funded  by  the  Naval  Material  Command  (08T23) , 

Task  Area  ZF43-421-001,  Work  Unit  1-1500-104-69. 

INTRODUCTION 

This  report  describes  and  presents  the  results  of  an  experiment  designed  to 

determine  the  unsteady  lift  and  drag  response  of  a  hydrofoil  running  in  regular 

waves  at  low  chord  Froude  numbers.  It  is  believed  that  these  experiments  are 

unique  in  that  they  were  run  at  low  Froude  numbers,  at  the  corresponding  frequency 

range  of  interest,  using  a  relatively  large  (0.406  meter  chord  and  2.43  meter  span) 

hydrofoil  model  and  in  that  tare  data  was  removed  to  isolate  the  foil  only  response 

A  literature  search  revealed  no  other  low  Froude  number  test  data. 

The  measured  lift  response  is  compared  to  the  prediction  obtained  from 

several  theories.  Extensive  theoretical  work  concerning  the  prediction  of  unsteady 

loads  on  air  foils  is  available.  References  1,  2,  and  3  summarize  and  list  the  out 

standing  works  applicable  to  the  prediction  of  airfoil  response.  In  this  report 

4  5 

the  theoretical  predictions  of  Theodorsen  (2-dimensional) ,  Reisner  ,  Reisner  and 

6  7  8 

Stevens  (3-D)and  Lawrence  and  Gerber  (3-D)  as  transformed  by  Henri  and  Ali  from 

the  case  of  forced  oscillations  in  an  infinite  medium  to  the  response  in  regular 

waves,  are  compared  to  the  measured  data.  In  addition,  a  lifting  surface  theory 

9 

prediction  as  computer  programmed  by  Widnall  and  computed  by  Besch  and  Food  for  a 
foil  very  similar  to  the  one  tested  here  is  compared  to  the  experimental  lift 
response. 

^References  are  listed  on  page  24. 


The  U.S.  Navy's  interest  in  the  test  data  presented  here  arises  from  possible 
design  applications  for  control  foils  (pitch  and  roll)  on  SWATH  (Small  Waterplane 
Area  Twin  Hull  Ship)  and  on  conventional  surface  ships  operating  at  slow  speed. 

In  addition,  these  tests  may  be  of  use  in  the  design  of  pure  hydrofoil  ships  with 
very  large  chord  foils.  For  example  the  test  covers  the  full  scale  speed  range  from 
18.3  knots  to  55.3  knots  for  a  hydrofoil  having  a  20  ft  (6.1  m)  mean  foil  chord. 

OBJECTIVES 

The  object  of  these  experiments  was  to  determine  the  lift  and  drag  response  of 
a  hydrofoil  operating  in  regular  waves  at  low  chord  Froude  numbers,  and  to  compare 
the  experimental  response  with  the  theoretical  response  using  existing  gust  response 
theories  available  from  aeronautical  or  marine  literature.  The  development  of  any 
new  response  theory  was  not  intended  even  though  it  was  recognized  that  the  available 
prediction  methods  were  not  likely  to  take  into  account  the  low  Froude  number  effects. 

DESCRIPTION  OF  EXPERIMENT 

The  experiments  were  conducted  in  the  Deep  Water  Basin  of  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center.  The  Deep  Water  Basin  is  2775  ft. 

(845.8m)  long  51  feet  (15.5  m)  wide  and  22  ft.  (6.70  m)  deep.  The  model  was  attached 
to  the  vertical  rails  of  Towing  Carriage  2  and  the  Pneumatic  Wavemaker  was  used  to 
generate  regular  waves.  For  the  wave  portion  of  the  tests,  the  water  level  was 
lowered  to  a  depth  of  19.5  ft.  (5.94  m) .  Additional  information  on  the  towing 
facilities  is  contained  in  Ref.  10. 

The  foil  model  used  in  these  experiments  has  been  extensively  tested  in  calm 
water  to  determine  the  effect  of  low  chord  Froude  number  on  the  wave  making  drag. 

These  calm  water  test  results  are  expected  to  be  published  subsequent  to  this  report. 
Figure  1  shows  a  photograph  of  the  foil  and  Table  1  lists  the  principal  dimensions 
and  other  important  information  concerning  the  foil. 

The  general  layout  of  the  test  set-up,  Figure  1  shows  the  tilt  yaw  table  which 
can  move  on  the  vertical  rails  for  height  adjustment,  a  forward  strut,  aft  strut 
and  the  top  pipe  which  provide  a  structural  base  for  block  gauge  stacks,  forward 
and  aft  fairings  around  the  block  gauge  stacks,  and  lower  pipe  and  nose  piece 
combination  which  supports  the  foil.  The  fairings  were  removed  in  Figure  2  to  show 
details  of  the  forward  block  gauge  stack  and  the  aft  block  gauge  flexure  combination. 
In  Figure  2,  the  foil  has  been  removed  and  the  nose  piece  slot  filled  with  a  faired 
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wooden  block  for  the  purpose  of  making  tare  measurements.  In  addition,  a  faired 
tail-cone  fitted  when  making  test  runs  has  been  removed  from  the  aft  end  of  the 
sting.  This  tail  cone  is  approximately  1  ft.  (0.304  m)  long  and  matches  the  6  inch 
(15.24  cm)  outside  diameter  of  the  sting. 


TABLE  1 

Dimensions  and  Characteristics 


Foil  Principal 

Fo  il  span 
Foil  chord 
Planform  shape 
Planform  area,  A 

P 

Aspect  ratio  AR 

Foil  material 

Foil  incidence  angle  at 
zero  speed 

Flap  angle 

Section  shape 

Trip  wire 


8  ft.  (2.438  m) 

1.33  ft.  (0.406  m) 
Rectangular 
10.66  ft.2  (0.99  m2) 
6 

Solid  aluminum 


NACA  64  A010 

0.016  inch.  (0.04  cm) 


The  primary  measurements  taken  and  the  instruments  used  for  taking  these 
measurements  are  shown  in  Table  2. 


TABLE  2 


Primary  Test  Measurements 

Forward  lift  5000  lb  (22.24  kN)  capacity, 

DTNSRDC  4  inch  (10.2  cm)  block  gauge 

Aft  lift  2000  lb  (8.89  kN)  capacity, 

DTNSRDC  4  inch  (10.2  cm)  block  gauge 

Drag  200  lb  (889  N)  capacity, 

DTNSRDC  4  inch  (10.2  cm)  block  gauge 


Wave  height  Wesmar  sonic  probe 

Carriage  speed  DTNSRDC  magnetic  pickup 
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The  forward  lift  and  drag  block  gauge  stack  was  located  3  feet  (0.914  m)  behind 

the  foil  \  chord  position  and  the  aft  lift  gauge  was  9  ft.  (2.74  m)  aft  of  the 

\  chord  position.  A  plan-view  of  the  foil  sting  set  up  is  shown  in  Figure  3.  A 

total  of  four  basic  foil/sting  test  configurations  were  tested  as  shown  in  Table  3. 


TABLE  3 

Test  Configurations 


Configuration 

1)  sting  only 

2)  sting  only 

3)  sting  4-  foil 

4)  sting  4-  foil 


Type  of  run 

Calm  water  tare 

Regular  wave-tare  runs, 
head  seas 

Calm  water 

Regular  waves  -  head  seas 


Water  depth 
22  ft.  (6.71  m) 
19.5  ft.  (5.94  m) 

22  ft.  (6.71  m) 
19.5  ft.  (5.94  m) 


By  subtracting  the  tare  runs  from  the  sting  4-  foil  runs,  the  lift  and  drag  of 
the  foil  only  were  isolated.  In  order  to  run  in  waves,  the  water  depth  had  to  be 
lowered  from  the  calm  water  level.  The  difference  in  water  depth  is  not  expected  to 
introduce  any  measurable  error  into  the  tare  subtraction  process. 

The  matrix  of  test  runs  involved  a  variation  of  foil  depth,  foil  forward  speed, 
wave  period  and  wave  height.  Foil  depth  was  fixed  for  each  run  at  foil  depth  to 
foil  chord  ratios  of  0.5,  1.0  and  2.0.  Foil  speed  varied  from  8.0  to  24.0  fps 
(2.43  to  7.32  m/s)  corresponding  to  a  chord  Froude  numbers  of  1.22  to  3.66.  Wave 
period  varied  from  1.75  to  4.25  sec.  The  nominal  wave  amplitude  for  the  test  was 
3.5  inches  (8.89  cm)  with  linearity  checks  at  nominal  wave  amplitudes  of  2.0  inches 
(5.08  cm)  and  5.0  inches  (12.7  cm). 

The  primary  test  speeds  corresponded  to  the  chord  Froude  numbers  F^c  of  1.22, 
1.83,  2.60  and  3.66.  Additional  speeds  were  run  as  required  in  order  to  obtain 
test  data  where  the  reduced  frequency  of  encounter,  k^,  was  held  constant.  Dynamic 
lift  and  drag  data  for  k^  =  1.83,  2.50,  and  3.14  were  obtained  for  a  range  of  chord 
Froude  numbers. 


EXPERIMENTAL  RESULTS 


* 


The  force  response  in  waves  was  harmonically  analyzed  to  determine  the 
amplitude  and  phase  at  the  fundamental  and  twice  the  fundamental  frequency.  The 
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measured  lift  and  drag  response  was  expressed  in  the  following  form: 
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where 


CL1’  CL2 


is  the  total  measured  lift  coefficient 

is  the  mean  value  of  the  lift  response  coefficient 

is  the  amplitude  of  the  oscillatory  part  of  the  response  at 
the  fundamental  and  twice  the  fundamental  frequency 


$  1 ,  is  t^ie  phase  angle  of  the  oscillatory  part  of  the  lift  response 

at  the  fundamental  and  twice  the  fundamental  frequency 

Cp  is  the  total  measured  drag  coefficient 

C  is  the  mean  value  of  the  drag  response  coefficient 


"Dl ’  D2 


is  the  amplitude  of  the  oscillatory  part  of  the  response  at  the 
fundamental  and  twice  the  fundamental  frequency 


$Di »  $D2  is  the  Phase  angle  of  the  oscillatory  part  of  the  fundamental  and 
twice  the  fundamental  frequency  drag  response 


The  subscript  1  and  2  associated  with  the  force  coefficients  and  phase  angles 
indicate  the  value  at  the  fundamental  frequency  or  twice  the  fundamental  frequency. 
Assuming  that  is  zero,  the  measured  lift  with  respect  to  the  incoming  wave  then 

can  be  shown  as  in  Figure  4,  In  order  to  isolate  the  foil  only  response  from  the 
combination  foil  and  sting  dynamometer  response,  tare  values  (sting  only)  were 
subtracted  from  the  magnitude  and  phase  of  the  response  for  the  sting  and  foil 
combination  at  both  the  fundamental  frequency  and  twice  the  fundamental  frequency. 
For  example,  for  determining  the  fundamental  frequency  lift  response  of  the  foil 
only  configuration  the  following  equations  are  used. 
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The  subscripts  f,  s,  and  fs  are  defined  below 

s  -  measured  response  of  sting  only  (tare)  configuration 
fs  -  measured  response  of  configuration  with  sting  and  foil 

f  -  response  for  the  foil  only  condition  as  calculated  in  Equations 
(3)  and  (4) 

The  drag  response  for  the  foil  only  condition  at  the  fundamental  frequency  and 
twice  the  fundamental  frequency  was  calculated  using  Equations  3-5  by  substituting 
drag  coefficients  and  phases  for  the  lift  coefficients  and  phases.  The  subscript 
"f"  indicating  a  foil  only  condition  will  be  dropped  from  here  on  in  this  report 
and  it  will  be  understood  that  the  tare  values  have  been  removed  from  all  the 
response  data  presented.  A  tabulation  of  the  unsteady  experimental  data  corrected 
for  tare  values  is  presented  in  Appendix  A. 

LIFT  AND  DRAG  RESPONSE  AT  NOMINAL  WAVE  HEIGHT 

The  amplitude  and  phase  of  the  lift  response  at  the  fundamental  frequency  for 
three  different  depth  chord  ratios  are  shown  in  Figures  5A  through  5F.  Almost  all 
of  the  data  was  collected  at  a  nominal  wave  amplitude  to  semi-chord  ratio  r  equal 
0.44,  In  order  to  maintain  clarity.  Figure  5D  shows  phase  data  for  only  selected 
chord  Froude  numbers.  The  data  for  other  chord  Froude  numbers  is  tabulated  in 
Appendix  A.  The  few  dashed  symbols  indicate  that  the  data  was  collected  at  some 
other  wave  height  at  which  good  linearity  was  observed.  Figure  6A  through  6F  and 
7A  through  7F  show  the  drag  response  and  phase  at  the  fundamental  and  twice  the 
fundamental  frequency. 


A. 
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The  data  for  the  lift  response  at  twice  the  fundamental  frequency  is  not  shown 
because  the  amplitude  of  the  response  C ^  was  very  small.  Typically  was  0.01 
to  0.02  times  the  primary  frequency  response  C^.  Almost  all  of  the  energy 
dissipated  in  lift  response  went  into  the  fundamental  frequency  response.  Values 
of  R00TQ0  for  virtually  all  of  the  test  runs  ranged  from  0.96  to  1.00.  As  mentioned 
in  Appendix  B  of  this  report  describing  the  sample  output,  R00TQ0  is  equal  to  the 
standard  deviation  of  the  oscillating  part  of  the  response  multiplied  by  .  In 
our  case  a  R00TQ0  value  equal  to  1.00  indicates  that  the  response  is  a  perfect 
sine  wave  of  the  fundamental  frequency.  The  R00TQ0  value  can  be  considered  as  a 
measure  of  the  energy  in  the  response  concentrated  at  or  near  the  fundamental 
frequency,  and  ROOT QO  values  near  unity  indicate  that  the  energy  of  the  response 
is  concentrated  at  the  fundamental  frequency. 

In  contrast  to  the  nature  of  the  lift  response,  the  drag  response  showed  a 
strong  component  at  the  second  harmonic.  values  as  high  as  70  percent  of  the 

Ch,  values  were  recorded.  Compare  Figure  6A  and  7A  at  k  =0.5  and  F  =1.22. 

Dl  t-a  e  nc 

This  strong  second  harmonic  drag  component  should  not  be  surprising  since 
reference  11  states  that  under  certain  condit ions  the  amplitude  of  the  second 
harmonic  drag  response  may  be  50  percent  of  the  fundamental.  A  tabulation  of  the 
lift  and  drag  response  data  appears  in  Appendix  A  grouped  according  to  depth  chord 
ratio.  Tare  values  have  been  substracted  from  all  the  response  data  and  the  run 
numbers  correspond  to  the  foil  and  sting  dynamometer  runs. 

In  general  the  phase  response,  especially  for  the  first  and  second  harmonic  of 
the  drag,  has  more  scatter  than  the  amplitude  of  the  response.  One  should  keep 
in  mind  however,  that  where  the  amplitude  of  the  response  is  small  or  near  zero  the 
associated  phase  is  a  meaningless  quantity  that  may  be  random. 

LINEARITY  OF  THE  RESPONSE 

The  bulk  of  the  test  data  was  taken  at  a  nominal  wave  amplitude  to  semi-chord 
ratio,  r,  equal  to  0.44.  Linearity  checks  at  nominal  values  of  r  =  0.22  and 
r  =  0.6  were  run  at  three  chord  Froude  numbers  and  corresponding  constant  reduced 
frequency  of  encounter  k^.  At  all  three  depth  chord  ratios  the  linearity  of  the 
fundamental  lift  response  is  very  good  as  shown  by  the  almost  horizontal  lines  in 
Figures  8A  to  8C.  The  linearity  of  the  amplitude  is  better  than  that  of  the  phase. 
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The  drag  response,  in  contrast  to  the  lift  response,  shows  poor  linearity 
(see  Figures  9A  to  IOC).  One  should  keep  in  mind  that,  as  stated  in  the  accuracy 
section  of  this  report,  the  drag  response  measurement  is  much  less  accurate  than 
the  lift  response  measurement.  In  addition,  the  drag  response  is  much  smaller 
than  the  lift  response  and  therefore  the  linearity  is  more  susceptible  to  the 
influence  of  scatter. 

RESPONSE  TRENDS  WITH  CHORD  FROUDE  NUMBER,  DEPTH  FROUDE  NUMBER  AND  DEPTH-CHORD  RATIO 

a)  Chord  Froude  number  (F  )  variation  at  constant  reduced  frequency  of 
encounter  and  depth. 

The  lift  and  drag  response  obtained  from  the  faired  lines  through  the  data 
shown  in  Figures  5A,  5C,  and  5D  were  plotted  holding  kg  and  depth  as  a  constant 
parameter.  See  Figures  11A,  11B,  and  11C.  The  constant  k^  values  0.183,  0.25  and 
0.315  were  selected  prior  to  the  experiment  and  the  proliferation  of  data  at  these 
encounter  frequencies  was  planned.  One  can  see  from  Figures  11A,  11B,  and  11C 
that  the  amplitude  of  fundamental  frequency  lift  response  strongly  increases  as 
chord  Froude  number  Fnc  decreases  from  3.66  to  1.22.  This  trend  agrees  with  other 
experimental  lift  response  trends  taken  at  high  chord  Froude  numbers.  (The  AG  (EH) 
data  shows  a  C^/r  change  from  0.065  to  0.092  for  a  chord  Froude  number  change  from 
7.18  to  A. 98,  see  Figures  2b  and  27  of  this  report).  The  data  for  the  fundamental 
drag  response.  Figures  12A,  12B,  and  12C  show  a  different  trend  with  chord  Froude 
number  variation.  The  fundamental  frequency  drag  response  curves  show  a  general 
hollow  at  Fnc  =  2.5  instead  of  the  steadily  decreasing  nature  of  the  lift  response 
curves.  The  second  harmonic  drag  response,  Figures  13A,  13B,  13C,  shows  the 
same  trend  as  the  fundamental  lift  response,  that  is,  decreasing  response  amplitude 
with  increasing  chord  Froude  number.  While  it  would  be  of  academic  interest  to 
find  an  explanation  for  this  trend,  in  any  practical  application  the  variation  in 
the  drag  response  would  probably  be  ignored,  since  its  value  is  very  small. 

The  solid  point  on  Figure  11B  represents  a  theoretical  prediction  for  a 
similar  foil.  This  prediction  is  discussed  in  the  theory  section  of  this  report. 
The  two  data  points  on  Figure  11C  were  taken  from  Reference  8  and  represent  results 
on  an  aspect  ratio  4  rectangular  hydrofoil.  They  are  shown  for  comparison  with 
our  aspect  ratio  6  hydrofoil  lift  response  trend. 
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b)  Chord  Froude  number,  F  ,  variation  at  constant  depth  Froude  number,  F  ,  , 

nc  nh 

and  reduced  frequency  of  encounter  k^. 

Figures  14A,  14b,  and  14C  shows  the  variation  of  the  fundamental  lift  response 
with  chord  Froude  number  holding  depth  Froude  number  and  constant.  These  curves 
are  all  in  a  narrow  band  indicating  that  the  response  is  sensitive  to  chord  Froude 
number  and  not  so  sensitive  to  depth  Froude  number.  For  example,  at  k^  =  0.25  and 
F  =  2,4,  Figure  14B,the  change  in  C  /r  corresponding  to  the  F  change  from  1.72 

HC  Li  TIC 

to  2.58  is  only  0.22  to  0.26.  Furthermore,  where  the  curves  cross  there  is  no 
dependence  on  F.  . 


c) 


Variation  of  the  fundamental  lift  response  with  F  ^  at  constant 
and  h/c . 


F 

nc 


A  carpet  plot  of  the  fundamental  lift  response  vs  F  ,  with  F  and  h/c  held 

nh  nc 

constant  is  shown  in  Figures  15A,  15b  and  15C  for  three  different  reduced  frequencies 
of  encounter.  The  symbols  shown  are  data  points.  These  plots  show  a  strong 
dependence  of  the  lift  response  on  chord  Froude  number  and  depth.  Note  that  as  k^ 
decreases  the  sensitivity  of  the  response  to  depth  Froude  number  diminishes. 


MEAN  LIFT  AND  DRAG  IN  WAVES  AND  IN  CALM  WATER 

The  mean  lift  and  drag  coefficients  C  and  C  in  head  seas  is  presented  in 

Li  1/ 

Figures  16A  to  16G  along  with  the  lift  and  drag  coefficient  and  in  smooth 
water.  All  these  coefficients  include  turbulence  stimulator  drag.  For  almost 
all  of  the  depth  chord  ratios  and  chord  Froude  numbers  the  mean  drag  in  waves  is 
less  than  the  smooth  water  drag  while  the  added  lift  in  waves  is  dependent  on  the 
depth-chord  ratio,  reduced  frequency  of  encounter  and  the  chord  Froude  number. 

In  order  to  assess  the  significance  of  the  lift  and  drag  response  in  terms 
of  the  mean  lift  and  drag  in  waves  the  ratio  and  are  plotted  in 

Figures  17Ato  17F.  These  figures  show  that  the  amplitude  of  the  fundamental  lift 
response  as  compared  to  the  mean  lift  is  always  much  larger  than  the  amplitude  of 
the  drag  response  compared  to  the  mean  drag.  The  highest  C^/C^  ratio  is  2.25  and 
occurs  at  depth  chord  ratio,  0.5,  reduced  frequency  of  encounter,  0.43,  and  chord 
Froude  number,  1.22.  At  these  same  conditions  the  C^/Cp  ratio  is  only  1.1.  In 
fact  the  amplitude  of  the  fundamental  frequency  drag  response  is  very  small,  probably 
negligible  for  most  control  surface  applications  at  depth  chord  ratio  1.0  and  2.0, 
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but  in  contrast  the  amplitude  of  the  lift  response  is  still  significant,  1.0  to 

0.8  times  the  mean  lift  for  F  =  1.22. 

nc 

THEORETICAL  PREDICTION  OF  LIFT  RESPONSE 

The  aerodynamic  response  of  wings  has  been  under  investigation  for  many  years. 

Several  texts,  (References  1,  2,  and  3)  deal  with  the  more  encompassing  aeroelastic 

problems  encountered  on  wings  and  they  should  be  consulted  for  background  information 

and  extensive  bibliographies.  There  appear  to  be  two  general  approaches  to  predicting 

the  slightly  different  problem  of  hydrofoil  response  in  waves.  One  approach  relies 

on  the  existing  aerodynamic  theories  (and  supporting  experimental  data)  as 

discussed  in  the  above  texts  for  predicting  the  aerodynamic  response  of  a  wing  due 

to  forced  sinusoidal  oscillations  in  infinite  uniform  flow  and  transforms  this 

response  to  that  of  a  hydrofoil  operating  under  regular  waves.  Examples  of  this 

A  5 

approach  include  response  predictions  by  Theodorsen  ,  or  Reisner  or  Lawrence  and 
Gerber^  transformed  to  the  hydrofoil  response  in  waves  by  the  approximate  method  of 

g 

Henry  and  Ali  .  The  other  general  approach  is  a  direct  attack  on  the  problem  of  a 

12 

hydrofoil  moving  under  regular  waves.  These  methods  such  as  Crimmi's  and 
13 

Widnalls  2-D  theory,  are  potential  flow  problems  which  satisfy  the  Laplace  equation, 

associated  kinematic  and  dynamic  free  surface  boundary  conditions,  and  the  Kutta 

condition  for  finite  pressure  at  the  trailing  edge  of  the  foil.  A  lifting  surface 
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program  developed  by  Ashley  et  al  and  Widnall  and  computer  programmed  by  Widnall 

9 

is  described,  summarized  and  evaluated  by  Besch  and  Rood  .  Further  details  and  a 
comparison  of  predictions  with  the  experimental  results  reported  herein  follow. 

a.  THEODORSEN' S  THEORY  (References  4,  and  16) 

This  is  a  2-dimensional  inviscid,  incompressible  linearized  potential 
flow  solution  satisfying  the  Kutta  condition.  The  foil  of  semi  chord  b  is  assumed 
to  be  in  uniform  flow  traveling  at  a  constant  speed  U  and  oscillating  with  a 
circular  frequency  w.  The  prime  motivation  for  the  theory  was  to  predict  flutter 
speed  and  the  lift  and  moment  response  on  airplane  wings. 

Results  are  expressed  as  a  tabulation  of  the  reduced  frequency  k  =  wb/V  versus 
a  function  C(k)  =  F  +  iG.  Since  the  publication  of  the  original  paper,  the  function 
C(k)  has  become  known  as  Theodorsen's  function. 

It  is  convenient  to  express  the  uniform  flow  response  to  forced  oscillations 
in  terms  of  the  non-dimensional  aerodynamic  coefficients  L^  and  L^.  They  are 
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related  to  Theodor sen's  function  by  the  following  expressions 


Lh  *  1  -  T  (m0> 


La  =  j  -  |  [1  +  2  (F+iG)]  -  -2  (F+iC) 

k 


and  L&  are  then  separated  into  their  real 
and  imaginary  parts: 


=  L, _ +  iLv  .  and 


"hr  hi 


L  =  L  +  iL  . 

a  ar  ai 


(6) 

(7) 

(8) 
(9) 


The  Theodorsen  2 -dimensional  prediction  of  these  aerodynamic  coefficients 
as  a  function  of  k  is  presented  in  Figures  18A  to  18D.  It  should  be  noted  here 
that  the  Theodorsen  function  C(k)  is  equal  to  an  algebraic  combination  of  the 
Bessel  function  Jq,  and  Y  ,  and  could  be  calculated  rapidly  by  computer. 
However,  since  Reference  6  provides  a  very  extensive  tabulation  of  Theodorsen' s 
function  it  was  not  necessary  to  write  a  computer  program  for  performing  this 
calculation. 

For  the  case  of  the  hydrofoil  oscillating  in  simple  harmonic  motion  defined  by 

,  iwt 

h  =  h  e  and 

o 

i(wt+3>  ) 
a  =  a  e  1 

o 

the  expression  for  lift  forces  per  unit  span  is  given  by 


(10) 

(11) 


L'  =  ffpb3  u>2  [Lh  £  +  a(La  -(h  +  a))Lh]  (12) 

where  h,  a  and  a  are  defined  in  Figure  19  and  represents  the  phase  lag  of 
the  pitch  motion  relative  to  the  heave  motion. 

Equation  12  predicts  the  response  for  an  oscillating  foil  in  calm  water.  It 
was  not  used  to  compare  our  experimental  results  with  the  theoretical  predictions, 


because  all  of  the  dynamic  experiments  were  run  in  waves  with  the  foil  fixed 
(within  the  limits  of  the  dynamometer  stiffness).  The  equation  is  presented  in 
order  to  illustrate  that  the  total  lift  response  has  contributions  from  both  the 
and  aerodynamic  derivatives.  The  comparison  of  these  functions,  as  predicted 
by  this  and  other  theories,  is  shown  in  Figures  18A  to  18D. 

In  order  to  predict  the  lift  response  for  the  fixed  hydrofoil  operating  in 

g 

regular  waves  the  theory  of  Henry  &  All  is  used  instead  of  Equation  12.  This 
theory  obtains  the  approximate  gust  response  operator  R  (kg)  from  the  non-dimensional 
aerodynamic  coefficients  and  L^.  It  assumes  a  fluid  vertical  velocity  distribu¬ 
tion  given  by  the  following  equation: 

v(x,y,t)  =  (y)rui  exp(-  sin  ^y)[x  -  (V  +  c)  t]  (13) 


It  shows  that  by  the  application  of  certain  simple  transformations  the  vertical 
velocity  distribution  on  a  fixed  airfoil  in  a  sinusoidal  gust  can  be  made  nearly 
the  same  as  that  on  an  oscillating  airfoil  in  uniform  flow.  Figure  20  shows  the 
definition  of  the  terms  used  in  Equation  13. 

The  wave  elevation  is  given  by: 


C(x,t)  =  (c/2)r  ccs  [x  -  (V+c)  t] 


(14) 


The  wave  celerity  c,  dimensionless  wave  number  K  and  wave  length  X  are  given  by 

c  =  g/w  (15) 


K  =  ttc  /  X 

X  =  (±n  deep  water) 


(16) 

(17) 


Henry  and  All's  theory  constrains  the  gust  response  operator  R(ke)  to  be  of 
the  form  given  by 


R(k  )  cos  m  t 
e  e 


UL1 

—  cos  (w  t  +  <J>  ) 

r  e  LI 


(18) 
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The  response  to  a  sinusoidal  gust  is  sinusoidal  at  the  same  fundamental  frequency 
with  a  phase  shift  4>  .  No  harmonics  are  permitted.  C^,  the  amplitude  of  the 

oscillating  part  of  the  lift  response  is  defined  previously  in  Figure  6,  and  r  is 
the  wave  amplitude  (1/2  the  crest  to  trough  distance)  to  semi  chord  ratio 


r  = 


(19) 


The  predicted  response  in  terms  of  the  aerodynamic  coefficients  is 


JL1  ,  ,  -kd 

-  =  uk  k  e 

r  e 


(X"  +  Y  ) 


(20) 


d  =  2h/c 


(21) 


where 


X  = 


hr 


(1  +  7  )  +  K  L  , 
k  ai 

e 


(22) 


Y  = 


L  -  (1  +  £  )  L.  , 
ar  k  hi 

e 


(23) 


and  the  phase  of  the  response  is 


$ 


LI 


tan 


(Y/X) . 


(24) 


The  reduced  frequency  of  encounter  k^  for  head  waves  is  defined  as 


k 

e 


,  uV 

2V  U 


(25) 


with 


(26) 


13 


The  phase  angle  is  the  phase  lead  of  the  lift  time  history  (at  midchord) 
with  respect  to  the  wave  time  history. 

It  is  important  to  note  that  when  the  foil  is  oscillating  in  uniform 
flow  the  forcing  frequency  is  characterized  by  the  reduced  frequency  k  whereas 
when  it  is  operating  in  waves  the  forcing  frequency  is  characterized  by  the 
reduced  frequency  of  encounter,  k^,  in  order  to  take  into  account  the  effect  of 
wave  celerity. 

b.  REISNER'S  THEORY  (References  5  and  6) 

Reisner  developed  a  theory  to  take  into  account  the  effect  of  finite  span 
on  the  airload  distribution  of  an  oscillating  wing  of  aspect  ratio  >3.0.  The 
final  results  of  this  theory  are  formulas  for  the  spanwise  distribution  of  air 
forces  and  moments  on  a  wing  of  finite  span.  In  these  formulas  the  effect  of  the 
three-dimensionality  is  a  correction  term  a  to  the  basic  function  C(k)  of  the 
two-dimensional  Theodorsen  theory.  Tabulated  results  of  the  correction  factor 
a  are  presented  for  wings  of  elliptical  planform  and  wings  of  rectangular 
planform.  In  both  cases  the  aspect  ratios  are  3  and  6  and  a  range  of  values  of 
the  reduced  frequency  k  are  covered. 

The  aspect  ratio  6,  rectangular  planform  correction  factor,  a  was  used  to 
calculate  the  aerodynamic  coefficients  and  L^.  The  tabulated  correction 
factors,  a,  as  presented  by  Reisner  are  for  the  semi-span  locations  0.0,  0.4,  0.8 
and  1.0  and  they  were  graphically  integrated  as  defined  by  Equations  29  through 
32. 

At  each  semispan  location 

Lh  =  1  '  T  (C(k)  +  0)  (29) 

K  =  h  (c  (k)  +  a)  -  -T  Ps  +  c<k>  +  a]  (30) 
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i  i . 


and  the  aerodynamic  coefficients  are 


b/2 


Lh  =  - 


4  * 


-b/2 


b/2 


L  =  ^  f  L' 

a  b  ]  a 


dy 


-b/2 


(31) 


(32) 


Symmetry  about  the  central  zero  semispan  station  is  assumed.  The  calculated 
coefficients  are  shown  in  Figures  18A  to  18D.  Once  the  aerodynamic  coefficients 

g 

L,  . ,  L,  ,  L  .,  and  L  are  determined,  the  Henri  and  Ali  theory  was  used  to 
hi  hr  ai  ar 

predict  the  response  in  waves  just  as  in  the  case  of  the  Theodorsen  theory. 

The  Reisner  theory  is  compared  to  experimental  results  for  depth-to-chord 
ratio  2.0  in  Figures  21  and  22.  Note  the  general  overprediction  of  the  amplitude 
of  the  lift  response  by  11  to  16  percent.  This  agreement  is  remarkable.  The 
trends  in  the  lift  amplitude  data  are  accurately  predicted.  The  phase  response 
is  overpredicted  by  as  much  as  40  percent;  however,  the  predicted  trend  with 
chord  Froude  Number  is  similar  to  that  obtained  in  the  experiments  (Figure  22). 

c.  LAWRENCE  &  GERBER  THEORY  (Reference  7) 

Lawrence  and  Gerber  developed  a  low  aspect  ratio  (AR  _<  4)  potential  flow 
theory  for  calculating  the  forces  on  a  harmonically  oscillating  wing  in  inviscid 
incompressible  flow.  The  theory  is  limited  to  wings  with  straight  trailing 
edges.  The  starting  point  for  the  theory  is  an  exact  two  variable  integral 
equation  for  an  oscillating  lifting  surface  in  incompressible  flow,  from  which 
an  approximate  integral  equation  in  one  variable  for  unsteady  flow  is  derived.  The 
aerodynamic  lift  and  moment  coefficients  for  the  limiting  case  of  aspect  ratio 
approaching  zero  are  shown  to  agree  with  the  Jones'*^  theory  extended  to  unsteady 
forces.  As  aspect  ratio  approaches  infinity,  the  two  dimensional  result  is 
obtained  for  rectangular  wings;  and  as  the  reduced  frequency  approaches  zero  the 
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steady  state  case  of  Reference  18  is  reproduced.  In  Reference  7  the  aerodynamic 
coefficients  are  tabulated  for  wings  with  rectangular  and  triangular  planforms 
of  aspect  ratio  0.25,  0.50,  1.0,  2.0,  and  4.0  for  a  range  of  reduced  frequencies 
k  from  0.0  to  1.0. 

In  order  to  obtain  results  for  aspect  ratio  equal  to  6.0,  the  aerodynamic 
coefficients  were  interpolated  from  the  rectangular  planform  aspect  ratio  1.0, 

2.0  and  4.0  values  as  predicted  by  Lawrence  and  Gerber  and  from  the  infinite  aspect 
ratio  two-dimensional  Theodorsen  results.  The  interpolated  coefficients  for 
aspect  ratio  6  are  shown  in  Figures  18A  to  18D  along  with  the  Reisner  results. 

Again,  in  order  to  obtain  the  response  in  regular  waves,  the  theory  of  Henry 

g 

6  Ali  is  applied.  Comparison  with  the  experimental  data  are  shown  in  Figures  23 
and  24 . 

The  theories  mentioned  so  far,  Theodor sen ' s,  Reisner's  and  the  Lawrence  and 
Gerber  theory  use  the  response  from  infinite  fluid  forced  oscillation  and  therefore 
would  not  be  expected  to  predict  the  effect  of  the  free  surface  on  the  response. 

Any  effect  due  to  the  presence  of  the  free  surface  could  be  the  result  of  a 
hydraulic  jump  over  a  very  shallowly  submerged  foil  or  to  a  modification  of  the 
wave  pattern  by  the  response.  According  to  Reference  19  there  is  no  uniform  agree¬ 
ment  as  to  whether  these  free  surface  effects  should  be  characterized  by  chord 
Froude  number  or  depth  Froude  number.  The  Henri  and  Ali  transformation  of  the 
aerodynamic  coefficients  to  response  in  regular  waves  take  into  account  only  the 
effects  due  to  wave  celerity.  In  order  to  take  Froude  number  effects  into  account 
in  the  Theordorsen,  Reisner  or  Lawrence  and  Gerber  theories  one  would  have  to 
modify  the  aerodynamic  coefficients. 

It  can  be  seen  from  Figures  21  and  23  that  the  amplitude  of  lift  response 
prediction  by  the  Lawrence  and  Gerber  theory  is  in  closer  agreement  with  the 
experimental  data  than  the  Reisner  theory  predictions,  although  the  difference 
is  probably  not  significant.  The  Lawrence  &  Gerber  theory  prediction  is  generally 
6.0  to  10.0  percent  too  high  at  the  depth  chord  ratio  2,0. 

One  would  expect  the  best  agreement  with  experimentally  determined  lift 
response  at  the  deepest  submergence  for  the  Reisner  and  the  Lawrence  and  Gerber 
theories,  since  they  do  not  take  free  surface  effects  into  account.  The  Lawrence 
6  Gerber  prediction  at  shallower  depths  are  shown  in  Figures  25A  to  25F.  Only  the 
Lawrence  &  Gerber  prediction  is  shown,  since  the  Reisner  lift  amplitude  predictions 
were  just  slightly  higher  for  the  depth  chord  ratio  2.0  case  and  one  would  expect 
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this  trend  between  the  two  theories  to  hold  independent  of  depth.  The  predicted 
trend  of  the  response  is  similar  to  that  obtained  from  the  experiments .showing 
decreasing  lift  amplitude  with  increasing  submergence  for  h/c  greater  than  1.0. 

However,  the  decrease  in  response  as  submergence  appraoches  zero  is  not  predicted. 

Note  also  that  the  predicted  phase  of  the  lift  response.  Figures  25D  to  25F,  is  con¬ 
stant  with  depth  chord  ratio  whereas  the  experimental  data  is  not. 

d.  WIDNALL  &  LANDAHL  THEORY 

13 

Widnall  and  Landahl  derived  an  acceleration  potential  theory  for  the  lift  and 
moment  response  of  an  oscillating  three-dimensional  hydrofoil  at  infinite  chord 
Froude  number  and  a  two  dimensional  theory  at  finite  chord  Froude  number.  The  two- 
dimensional  theory  contains  the  effects  of  the  unsteady  wake,  surface  waves  generated 
by  the  motion  and  the  depth  of  the  foil  below  the  surface.  At  infinite  depth, 
this  theory  is  claimed  to  reproduce  the  Theodorsen  results. 

e.  WIDNALL  LIFTING  SURFACE  PROGRAM 

A  lifting  surface  theory  for  calculating  the  fluid  dynamic  loading  on  airfoils 

14 

and  hydrofoils  was  developed  by  Ashley  et  al  and  Widnall  and  was  subsequently 

9 

programmed  into  computer  code  by  Widnall.  Besch  and  Rood  performed  extensive 
calculations  and  comparison  to  existing  data  using  this  program  and  have  shown  that 
with  empirical  corrections,  lift  predictions  can  be  made  with  an  accuracy  of  about 
15  to  20  percent.  The  computer  program  is  capable  of  calculating  both  steady  and 
unsteady  loading  on  hydrofoils  (camber  taken  into  account,  thickness  effects 
ignored)  and  includes  finite  depth  effects  and  cavitation  effects.  A  summary  of  the 
lifting  surface  theory,  a  program  listing  and  user's  guide  is  presented  by  Besch 
and  Rood^. 

9 

Besch  and  Rood  used  the  Widnall  program  to  predict  the  lift  response  for  an 
aspect  ratio  6  foil,  on  which  lift  response  tests  were  performed  by  Wetzel  and 
Maxwell^.  This  foil  was  similar  to  the  one  tested  here  except  that  it  has  NACA 
16-509  section  and  had  a  much  smaller  chord  =  0.25  ft  (7.62  cm).  It  has  the  same 
rectangular  planform  and  the  lift  response  prediction  was  made  for  kp  equal  to  0.237 
and  depth  chord  ratio  1.0  as  shown  in  Figure  12B.  The  prediction  for  the  Wetzell 
and  Maxwell  foil  which  is  very  similar  to  the  present  configuration  compares  well  j 

with  the  present  data. 

In  order  to  compare  the  various  prediction  methods,  the  Lawrence  and  Gerber  j 

and  the  Theodoresen  2-dimensional  theories  were  used  to  predict  the  lift  response  j  ] 
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amplitude  of  the  AG(EH)  main  foil.  Widnall  2-dimensional  and  3-dimensional 

predictions  as  well  as  experimental  data  are  available  for  this  foil  as  shown  in 

Figures  26  and  27.  The  AG(EH)  main  foil  is  a  trapezoidal  aspect  ratio  3  foil. 

Trapezoidal  planforms  can  be  predicted  by  the  Lawrence  and  Gerber  theory;  however, 

the  aerodynamic  coefficients  were  available  in  tabulated  form  only  for  the 

triangular  and  rectangular  planforms.  It  is  assumed  that  the  prediction  for  the 

trapezoidal  planform  would  fall  between  the  triangular  and  rectangular  planforms 

aspect  ratio  3  prediction.  It  can  be  seen  from  Figure  26  and  27  that  the  Widnall 

9 

3-dimensional  prediction  as  calculated  by  Besch  and  Rood  is  the  best.  As 
expected,  the  two-dimensional  Theodorsen  theory  greatly  overpredicts  the  response. 

The  Lawrence  and  Gerber  theory  also  over predict  the  response  but  not  by  much. 

Based  on  previous  trends,  the  Lawrence  and  Gerber  theory  would  probably  be  more 
accurate  at  a  deeper  submergence.  The  Reisner  theory,  not  shown, would  be  expected 
to  be  very  close  to  the  Lawrence  and  Gerber  theory.  An  anomaly  in  the  prediction 
is  the  underpredict  ion  by  the  two-dimensional  Widnall  theory  obtained  from 
Reference  20.  One  would  expect  the  two-dimensional  theory  to  overpredict  the 
experimental  response  as  well  as  the  three-dimensional  theory. 

DYNAMOMETER  CALIBRATION  AND  ACCURACY 

The  set-up  and  calibration  for  the  experiment  consisted  of  calibrating  the 
individual  block  gauges,  and  re-calibrating  and  checking  out  the  assembled  gauges 
using  the  Carriage  V  heavy  bridge  as  a  calibration  stand.  Then,  without  disassembly, 
the  dynamometer  was  moved  to  the  Carriage  II  dry  dock  and  a  final  calibration  was 
performed. 

Based  on  past  experience,  an  interaction  between  the  lift  and  drag  measurement 
was  expected.  This  interaction  is  basically  caused  by  small  deflections  in  the 
sting  and  sting  supports  when  they  are  under  a  lift  load  and  by  possible  misalign¬ 
ment  of  the  flexure.  The  interaction  was  minimized  by  using  shims  to  adjust  the 
vertical  angle  of  the  flexure  within  the  limits  allowed  by  the  clearance  in  the 
flexure  bolt  holes.  The  angular  deflection  of  the  nose  piece  at  the  1/4  chord 
position  was  determined  by  applying  varying  lift  loads.  The  final  deflection 

calibration  performed  on  the  assembled  dynamometer  in  the  Carriage  II  drydock 
-4 

showed  3.74  x  10  degrees  of  deflection  per  pound  (4.48  N)  of  lift  applied  at  the 
1/4  chord  position  of  the  foil. 
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The  excess  interaction  drag  produced  by  various  lift  loads  was  found  to  be 
consistently  repetitive  and  was  calibrated.  The  calibration  curve  was  fitted  with 
the  following  polynomial  expression 

AD  =  -0.035  +  0.00146  L  +  1.9375  x  10~6  L2  (30) 

where  AD  is  the  excess  interaction  drag  in  pounds  (4.448N)  due  to  the  applied 
lift,  L  in  pounds  (4.418  N)-  AD  is  equal  to  the  channel  1  mean  DRAG  minus  the 
channel  12  DRAG  COR  mean  drag  as  shown  on  the  sample  outputs.  Appendix  E. 

From  the  equation  for  AD  it  can  be  seen  that  the  maximum  interaction  drag  is 
10.7  lb.  (47.59  N)  for  2000  lb  (8896  N)  applied  lift.  The  accuracies  of  the 
primary  measurements  used  in  this  experiment  are  shown  in  Table  4. 

TABLE  4 

Primary  Measurement  Accuracies 
Measurement  Accuracy 

Lift  +  5%  of  reading  for  lift  <200  lb. 

+  1.5%  of  reading  for  lift  >200  lb. 

Drag  +  1.5%  of  reading  for  drag  <50  lb. 

+  0.5%  of  reading  for  drag  >f>0  lb. 

Speed  +  0.01  knot 

Wave  height  +  2%  of  reading 

Wave  period  +  0.05  sec. 

The  lift  and  drag  accuracies  were  determined  by  the  calibration  of  the 
assembled  dynamometer.  The  drag  accuracy  includes  the  correction  for  interaction 
drag  caused  by  lift.  The  wave  height  and  speed  accuracies  are  based  on  previous 
calibration  and  experience  with  the  instruments  and  the  wave  period  accuracy  is 
based  on  the  scan  rate  of  the  data  collection  system. 

The  above  instrument  accuracies  were  used  to  calculate  the  error  range  of  the 
measured  lift  and  drag  coefficients,  lift  and  drag  response  and  reduced  frequency 
of  encounter  for  two  typical  data  sets  which  were  chosen  to  represent  the  best  and 
the  worst  expected  errors.  These  errors  are  shown  in  Table  5  both  as  an  absolute 
value  and  as  a  percentage  error. 


TABLE  5 

Accuracy  of  the  Reduced  Data 


Low  Speed 

High  Speed 

F  -  1.22 
nc 

F  =3.66 
nc 

Nominal 

Value 

Error 

(absolute) 

%  Error 

Nomi  > . ' 
Va.1 

Error 

(absolute) 

%  Error 

Lift  coefficient 

0.202 

+0.011 

+5.4 

0.266 

+^.004 

+1.5 

Drag  coefficient  Cp 

0.0234 

+0.004 

+1.7 

0.0247 

+0. 00014 

+5.7 

Lift  response  CTl/r 
(fundamental  freq) 

0.4204 

+0.027 

+6.4 

0.165 

+0.010 

+6.1 

Drag  response  Cp^/r 
(fundamental  freq) 

Reduced  frequency  of 
encounter  k 

e 

0.0066 

+0.0013 

+19.7 

0.0034 

+0.00079 

+23.2 

Low  wave  period  at 

t  =  2.0 
o 

0.466 

+0.017 

+3.6 

0.291 

+0.012 

+4.1 

High  wave  period 

at  x  =  4.0 
o 

0.182 

+0.003 

+1.6 

0.0947 

+0.0018 

±!  .9 

It  can  be  seen  that  the  lift  coefficient  C  ,  drag  coefficient  C  ,  lift  response 

Lj  u 

C^/r,  and  reduced  frequency  of  encounter,  k^,  have  low  error  values.  The  drag 
response,  C^/r,  on  the  other  hand  has  a  very  large  error  +19.7%  and  23.2%  for  the 
two  cases  considered.  This  is  primarily  the  result  of  trying  to  measure  a  very 
small  alternating  drag  with  a  block  gauge  which  must  have  the  capacity  to  measure 
at  least  an  order  of  magnitude  larger  steady  state  force.  It  should  be  emphasized 
that  the  error  limits  in  Table  6  are  strictly  instrument  errors  and  as  such  they 
do  not  account  for  scatter  in  the  test  data  which  is  due  to  a  variation  in  test 
conditions  or  the  non  repeatability  of  any  possible  flow  that  is  in  the  transition 
stage  between  the  laminar  and  turbulent  regime. 

In  addition  to  static  calibrations,  the  dynamic  response  of  the  dynamometer 
(with  the  foil  mounted)  was  explored.  A  small  mechanical  shaker  shown  in  Figure  28 
was  used  to  impart  dynamic  lift  and  drag  forces  to  the  dynamometer  at  various 
frequencies.  In  order  to  measure  the  lift  dynamic  response  a  small  A-frame  was 
constructed  which  held  the  shaker  above  and  out  of  the  water  at  the  1/4  chord  foil 
positions  while  the  foil  was  submerged  to  a  depth  chord  ratio  1/2.  The  drag 
dynamic  response  was  taken  with  the  foil  raised  out  of  the  water,  the  tail  cone 
removed  and  the  shaker  bolted  in  its  place  so  as  to  produce  sinusoidal  horizontal 
forces  in  line  with  the  center  of  the  sting. 

The  drag  amplification  factor  was  taken  in  air;  however, it  is  felt  that  the 
results  would  be  nearly  the  same  in  water  because  the  added  mass  associated  with 
the  drag  direction  is  expected  to  be  small.  The  resonant  frequency  of  the  drag 
response  in  water  was  double  checked  by  striking  the  dynamometer  with  a  rubber 
mallet  and  analyzing  the  strip  chart  record  of  the  drag  gauge  output  as  it  decayed 
over  approximately  10  cycles.  The  above  method  yielded  a  4,8  hertz  natural 
frequency  in  water  versus  the  4.5  measured  in  the  air  shake  test.  Although  one 
would  expect  a  slightly  lower  natural  frequency  in  water  it  is  felt  that  for  the 
purposes  of  this  experiment  the  drag  amplification  factor  could  be  considered  the 
same  in  water  and  in  air.  The  results  shown  in  Figures  29  and  30  indicate  that  the 
amplification  factor  is  unity  for  both  lift  and  drag,  in  the  low  frequency  range, 
where  the  experiments  were  conducted.  Note  that  in  Figure  29  the  shaker  lift  force 
is  double  that  of  the  force  in  block  gauge  #2  for  the  low,  less  than  2.0  hertz, 
frequencies  because  by  static  analysis,  the  112  block  gauge  records  only  1/2  the 
lift  force.  The  maximum  encounter  frequency  run  during  the  tests  was  1.65  hertz 
while  the  bulk  of  the  test  data,  over  80%  of  the  runs,  were  at  encounter  frequencies 
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between  033  and  1.1  hertz.  Thus,  both  lift  and  drag  amplifications,  in  the 
frequency  range  of  interest  are  considered  to  be  1.0. 

The  lift  amplification  curve.  Figure  29,  has  a  double  peak  which  is  typical 
of  a  "vibration  absorber"  (see  Reference  21)  that  is  made  of  a  pair  of  idealized 
spring-mass  dampers  in  series.  The  analogy  could  be  made  to  the  test  set  up,  where 
the  struts  and  the  frame  supports  consist  of  one  spring  mass  system  and  the  foil¬ 
sting-block  gauge-combination  is  the  other  spring  mass  system. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Based  upon  experimental  results  and  theoretical  predictions  the  following 
conclusions  and  recommendations  can  be  made: 

CONCLUSIONS 

1.  Chord  Froude  number  has  a  strong  influence  on  the  amplitude  and  phase  of  the 
first  harmonic  lift,  first  harmonic  drag  and  second  harmonic  drag  response. 
Chord  Froude  number  effects  become  more  important  as  depth  chord  ratio 
decreases.  At  constant  reduced  frequency  of  encounter,  the  fundamental 
frequency  lift  response  increases  with  decreasing  chord  Froude  number  for  all 
depth  to  chord  ratios  and  chord  Froude  numbers  from  1.22  to  3.66  that  were 
tested. 

2.  The  second  harmonic  lift  response  is  zero. 

3.  The  amplitude  and  phase  of  the  lift  and  drag  response  vary  with  the  reduced 
frequency  of  encounter. 

A.  Of  the  three  depth  chord  ratios  tested,  0.5,  1.0  and  2.0,  the  amplitudes  of  the 
lift  and  drag  response  were  least  for  depth  chord  ratio  2.0. 

5.  The  first  harmonic  lift  response  amplitude  and  phase  are  linear  with  wave 
amplitude. 

6.  The  drag-response  amplitude  (both  first  and  second  harmonic),  in  general,  does 
not  appear  to  be  linear  with  wave  amplitude.  The  magnitude  of  the  measured 
oscillating  drag  component  was  very  small,  posing  difficulties  in  data- 
resolution.  Guidance  from  theoretical  predictions  is  needed  for  establishing 
data  trends. 

7.  Prediction  of  the  fundamental  lift  response  amplitude  at  depth  chord  ratio  2.0 
by  the  theory  of  Reisner  is  generally  within  11  to  16  percent  of  the 


experimental  data.  The  predictions  of  Lawrence  and  Gerber  for  fundamental  frequency 
lift  amplitude  are  within  6  to  10  percent.  These  predictions  do  not  take  into 
account  the  presence  of  the  free  surface.  As  depth-chord  ratio  decreases,  they 
deviate  more  and  more  from  the  experimental  data. 

RECOMMENDATIONS 

1.  Exercise  the  Widnall  lifting  surface  program  and  compare  the  lift  response 

predictions  to  the  present  experimental  data.  Concentrate  on  the  low  chord 

Froude  number  F  =  1.22,  low  depth  chord  ratio,  h/c  =  0.5  case, 
nc 

2.  Repeat  some  of  the  present  tests  with  a  geosim  model  of  smaller  size  (represent¬ 
ative  of  typical  SWATH  model  appendages)  in  order  to  determine  any  scale 
effects,  and  verify  Froude  number  effects  by  varying  the  chord. 

3.  Develop  a  drag  response  prediction  method  for  the  three  dimensional  case  that 
meets  the  need  for  accurate  prediction  at  low  chord  Froude  numbers  near  1.0  and 
at  low  depth  chord  ratios  near  0.5. 
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Figure  1  -  Photograph  of  Hydrofoil  &  Dynomometer  Arrangement 
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Figure  2  -  Details  of  Forward  &  Aft  Block  Gauge  Stacks 
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Figure  5D  -  Lift  Phase  Response  at  the  Fundamental 
Frequency;  h/c  =  1.0 
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Figure  5F  -  Lift  Phase  Response  at  the  Fundamental 
Frequency;  h/c  =  2.0 
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Figure  6A  -  Drag  Amplitude  Response  at  the  Fundamental 
Frequency;  h/c  =  0.5 
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Figure  6B  -  Drag  Phase  Response  at  the  Fundamental 
Frequency;  h/c  =  0.5 
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Figure  7D  -  Second  Harmonic  Drag  Phase 
Response;  h/c  =  1.0 


45 


h/c  =  2 


Figure  7F  -  Second  Harmonic  Drag  Phase 
Response;  h/c  =  2.0 
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Fundamental  Frequency  Lift  Response 
Linearity;  h/c  =0.5 
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Figure  8B  -  Fundamental  Frequency  Lift  Response 
Linearity;  h/c  =  1.0 
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Figure  9C  -  Fundamental  Frequency  Drag  Response 
Linearity;  h/c  =  2.0 
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Figure  14B  -  CT1/r  vs  F  ;  F  =  Constant;  k  =  0.250 
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Figure  14C  -  C  ,/r  vs  F  ;  F  =  Constant;  k  =  0.315 
Ll  nc  nn  e 
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Lift;  h/c  =  0.5;  F  =  1.22,  1.83  and  3.66 
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Figure  16D  -  Mean  Lift  in  Waves  Compared  to  Calm  Water 
Lift;  h/c  *=  1.0 


ke 

Figure  17B  -  Fundamental  Frequency  Drag  Amplitude  Response 
to  Mean  Drag  Ratio;  h/c  =0.5 
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Figure  17D  -  Fundamental  Frequency  Drag  Amplitude  Response 
to  Mean  Drag  Ratio;  h/c  =  1.0 
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Figure  17E  -  Fundamental  Frequency  Lift  Amplitude  Response 
to  Mean  Lift  Ratio;  h/c  =  2.0 
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Figure  18D  -  Theoretical  Prediction  of  k  L^r 
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Figure  22  -  Comparison  of  $L1  as  Predicted  by  Reisner's 
Theory  to  the  Present  Experimental  Results;  h/c  =  2.0 
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Figure  26  -  Comparison  of  Several  Theoretical  Predictions 
of  the  AG (EH)  Main  Foil  Lift  Amplitude  Response 
with  Experimental  Data;  F  =  4.98;  h/c  =  1.15 
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Figure  29  -  Lift  Response  Amplification  Factor 
(in  water,  h/c  =  0.5) 
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APPENDIX  A 

TABULATED  EXPERIMENTAL  DATA 

Tabulation  of  the  experimentally  determined  lift  and  drag  response  on  the 
aspect  ratio  6  hydrofoil  at  depth  chord  ratio  0.5,  l.C  and  2.0 

The  lift  and  drag  response  listed  in  this  appendix  corresponds  to  the  foil 
only  condition.  Dynamometer  tares  have  been  subtracted  from  the  initial  data. 


J 

* 


0.249  1.220  0.863  2.988  0.449  0.432  76.4  0.0213  1  71  0.00812 


301  0.296  1.837  1.300  I  5.345  0.271  0.349  71.8  0.0022  -59  0.00334 


333  0.464  1.217  0.861  5.515  0.249  0.539  78.7  0.0185  87  0.00870  -31 


317  0.171  3.040  3.040  5.116  0.445  0.191  68.4  0.0009  88  0.00070  -17 


APPENDIX  B 


SAMPLE  EXPERIMENTAL  OUTPUT 
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SAMPLE  EXPERIMENTAL  OUTPUT 

A  sample  of  the  data  output  obtained  for  each  test  spot  is  shown  in  Figure  IB 
and  2B.  The  output  consists  of  two  pages.  Page  1  lists  steady-state  type  and  mean 
value  information  and  page  2  lists  the  harmonic  analysis  for  unsteady  measurements 
recorded  in  wave  tests.  Page  2  is  deleted  during  calm  water  tests. 

The  experimental  conditions,  encounter  frequencies,  wave  lengths  and  water 
properties  are  shown  as  header  type  information  and  are  repeated  on  page  2  of  the 
sample  output. 

The  primary  measurements,  drag,  lift  forward,  lift  aft,  wave  height  (single 
amplitude),  and  carriage  speed  are  shown  as  raw  data  in  engineering  units,  on 
channels  1  through  6  of  the  output. 

Channels  7  through  12  represent  lift  and  drag  corrections  to  the  raw  data 
which  account  for  dynamometer  interaction  and  unwanted  lift  effects  due  to  slight 
pitch  angle  deflections  of  the  sting  dynamometer.  The  corrected  lift  and  drag 
are  shown  in  channels  10  and  12.  Note  that  channels  5,  8,  and  11  are  not  used. 

LIFT  AND  DRAG  CORRECTIONS  SHOWN  ON  SAMPLE  OUTPUT 

In  the  ideal  case  with  a  perfectly  rigid  dynamometer  there  would  be  no  need  to 
correct  the  measured  lift  and  drag.  However,  the  sting  support  deflects  slightly 
under  the  application  of  lift.  As  much  as  2000  lbs  (8896  m)  lift  was  generated 
during  some  of  the  runs.  The  slight  bending  of  the  sting  caused  two  problems: 

1)  A  repeatable  non-linear  dynamometer  interaction.  The  application  of  pure 
lift  results  in  a  small  drag  force. 

2)  A  discrepancy  between  the  totally  rigid  (no  angular  or  other  deflections) 
mode  of  operation  of  the  foil  as  assumed  by  theory  and  the  slight  elasticity 
in  the  foil  -  dynamometer  arrangement  of  the  test  set  up.  As  a  result  of 
large  lift  forces,  the  angle  of  attack  of  the  foil  deflected  from  its  4° 
static  incidence  angle  to  some  higher  mean  angle  of  attack.  During  foil 
operation  in  waves  the  angle  of  attack  oscillates  slightly  in  response  to 
the  varying  lift. 

The  problem  of  the  non  linear  interaction  drag  produced  *by  lift  was  solved  by 
static  calibration  as  discussed  under  the  calibration  section  of  this  report.  The 
amount  of  interaction  drag  is  the  difference  between  the  channel  1  drag  (which  is 
the  cal  factor  times  the  voltage  difference  of  the  drag  gauge)  and  the  channel 
12  DRAG  COR  value.  The  channel  12  DRAG  COR  is  the  value  of  the  drag  corrected  for 
dynamometer  interaction.  In  the  sample  case  of  Figure  IB  the  MEAN  drag  was  adjusted 


118 


CARR IAGc  II  LOU.  FPOUDE  NUMBER  CONTROL  FOIL  27  AUGUST  1973 


jL 


u  n 

a  -i 

Cu  I 
Cu 
•XCu 
-  J  i_ 
•X 


Mi 


119 


TiFf'TH  FHTiu  :  l.nnnn 


i_T| 

UJ 

LU 

*“■ 

to 

o'1 

— 

F 

LX' 

G'  CO 

i.n 

•— 

IX' 

rj  ~ 

Cl 

CO 

TO 

J  U"' 

K3 

TO 

•X 


I  - 
(M 


LL) 

F 


II  II 

>X 
F  \ 
<X 

L LM 
UJ  CD 
JZ 
LUX 
U  _J 


TO 
CD  1 


X 

X 

Cl 


Cvj 

Cl 


cri  K'  r  -  u'*  co  co  J  -  r  -  to  ud 

QU  CO  CTi  '-£>  Lj~i  U~  1  X  "T  *-* 
— ’  U~)  f  -  U"l  Ijj  *—  Lf»  'X  — 

I  —  '“•I 

I  I 

TO  TO  Cl  r  j  O  d  r-  j  O  TO  C 
O  Cl  C  C  C1  O  TO  TO  Cl  o 
+++1++I+++ 
UJ  UJ  UJ  UJ  LU  LlJ  iU  UJ  LU  LU 

vc  i  -  —  i  l*  '  'jt  co  u'i  at  c 

TO  *7  CO  K>  T  O  'X*  Oi  K» 

—  (‘  I  —  ••  •  CO  —  CO  t  _  TO  ~ 


CD 

¥ 

UJ 

<M 

c 


in 

OJ 

QD 

S3 

N 

A 

Lj 


Cl  ‘X  •— 
Ca  -7i  — 
(_i  ij,  IT ' 
C  TO 


CL  'XI  CO 
_J  X  IT) 
O  CJ  — 

co  r*  j 

ro  o 


«  f'J  i  -  TO 


r  —  X  t  • 


Cl 

t— i 
LJ 


l  M  C 
T(  U'i 


UT  C 
_J  <T,  K 
TO  'Li  f  - 


*-  X  .  T 
i“  —  X  _ 
L, _ ‘_ 


to  U-  Cl 


lD 

UJ 


J  X 
.  l  - 


J  X  <JJ 


V  * 

u  *  <x  Li"1  u  .  TO  IT ' 

WO 

to 

1 

lo  l  ' 

TOl  1 

to 

TO 

rsj 

00 

11 

TO 

m 

r  j 

(X  C-  J  T  - 

r-'i 

TO 

o 

Cl  TO  'J, 

o  TO  a, 

J  • 

■'X  -  - 

1  o  -  o  o  — 

O  C1 

TO 

C* 

i_i  '/  TO 

'•  J 

+  ^+ 

+  +  ++■  +  T 

+  + 

M 

ll 

. — 

LL!  X  LU  LU  LU  1X1  lx  Uj 

UJ  LU 

L;_ 

X 

—  lj. 
K>  ''-j 

|"  J  J-'  'Ll  ro  u  ‘  C 
(  -  U~|  J  CO  Li"  — 
J'j  t  u"'  -  r  - 

O  L 

J  ljj 
•■] 

FNH 

r 

LJ 

it- 

_J  r  -  r  _ 

r- 1  TO 

„  r  ,  CO  CO  -  O' 

CO  f  i 

Ci.  '-  '  5  f  - 
O  U~i  f  - 

')  TO  J  I 
4  4+’ 


.to  ix  j  uj 
“  r~" 


Uj  Uj  UJ  1 1  i  'i  J I J  ‘X  lx  lj  lx' 
■Xi  t  LT  Ci  In  r .  L"  TO  O', 
r.o  i  .  r  ,  »j.  [  _  '1  .o  r  -  ■  .  to 


C£| 

CT, 

K' 

T~t 

CT) 

Oj 

|‘VJ 

CO 

r*- 

• 

u 

Tj- 

TO 

TO 

TO 

II 

ii 

ii 

LJ 

_ J 

Lr. 

X 

•X 

LL 

•I 

i  -i.r,  T.  U .  .J .  >Xi  '.Ti.j.ijf  _ 

:r 

X 

LL 

1'"' 

»'  J  'X' 

L  ' 

U“ 

i  £i  r 

1 

to  —  “  ic  —  -.  i-: - -  o 

C  i—  x  x  i _ ■  ix  i  -  x  ix  TO 

TO 

4  rf  +  4  44+4  + 

•X 

'X 

U  '  *X» 

TO» 

LU  UJ  'U  LU  LJ  Lx  UJ  UJ  LU  LU 

c- 

h- 

'  7  L  i  "T  i  \J  f"  i 1  £i  »Ti  X  1 4 1  CO 

T  lx  TO  CO  *T  ,-'.j  'X  .  ■  J  ’ — < ' — ■ 

TO 

TO 

IX 

-- 

'X 

Ti  rx  'X  *X  F —  Lj  M  OOj 

" 

m  .■ 

LU 

:L 

lVj  -j-  —  —  — *  CO  — •  **•  j 

TO 

r- 

.X 

X  _5 

L 

ii 

« 

n 

ul  X 

'’•J 

LU 

■X 

LU 

LL 

, , 

- —  > — ■ 

LL 

IX  'XI 

f  i 

X  Lxi 

• — ■ 

i  i _ 

•  — 

IJ‘|  '  ■*'  I.  1  u  '  .  J  _l  1  I'l 

x  _ 

> - * 

LLI 

UjUULU  - - U  Tib  XI  _LX 

•I 

1 - 

X 

_J  — 1  -J  —  _l  —1  -I  —  _l 

*  r  - 

■  Tf 


LU  LU 
LX  LL 


.LL  # 


■X  LL  LL  ->  LL  —ILL  X 
Li.  X  —  —  LU  U. 


Figure  2B  -  Sample  Data  Output  -  Page  2 

120 


JHWJUIHP  R>  BjRT! 


from  the  16.62  lb  (73.92  N)  channel  1  value  to  the  corrected  value  16.45  lbs 
(73.16  N)  as  shown  In  channel  12.  Similarly  the  AMP1  and  AMP2  values  of  the 
oscillating  drag  were  adjusted  as  shown  in  Figure  2A. 

The  solution  to  the  second  problem,  that  of  the  foil  incidence  angle 
deflection  along  with  the  small  spanwise  bending  of  the  foil  and  the  small 
deflections  in  heave  constitute  a  hydroelastic  problem,  the  complete  solution  to 
which  is  outside  the  scope  of  this  experimental  effort.  However,  the  problem  was 
not  ignored.  It  is  assumed  that  the  foil  incidence  angle  deflection  will  have 
the  largest  effect  on  the  lift  response  and  a  quasi-steady  correction  method  was 
developed.  Details  of  the  method  follow. 

1)  Calculate  6a,  the  change  in  angle  of  attack  due  to  an  increment  in 
foil  lift. 

By  static  calibration: 

6a  =  3.74  x  l(f4  x  6L  (IB) 

where  6L  is  in  pounds  (4.448  N) 

6a  is  in  degrees 

Depending  on  the  harmonic  of  the  fundamental  frequency  under  consideration, 
6L  is  obtained  from  the  AMPl  or  AMP2  lift  response  value  found  in  channel  7. 


2) 


From  previous  unpublished  smooth  water  experiments  on  this  foil  obtain 
3C  /3a  and  3C  /3a  as  listed  in  Table  IB. 

Li  Is 


TABLE  IB 

2 

Variation  of  Lift  Curve  Slope  and  3CD/3a  with 
Depth  Chord  Ratio  and  Froude  Number 


!cl 

h/c 

3a 

3a2 

2.0 

0.038775 

(F  ) 
nc 

2  -  2.0347  (F  )  +  0.886 
nc 

0.00043 

1.0 

0.001718 

(F„r) 

nc 

+  0.0527 

0.00041 

0.5 

0.004206 

(F  J 
nc 

+  0.03682 

0.00038 
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3)  Calculate  AC^  and  AC^  -  the  lift  and  drag  coefficient  change  due  to  a 
small  change  in  angle  of  attack. 


AC 


L 


6a 


(2E) 


8CD  2 

ACd  =  — j  •  Oar 

3a 


(3B) 


4)  Calculate  AL  -  the  lift  and  drag  due  to  a  small  change  in  angle  of 
attack: 


AL  =  ACl  (h  PV2  Ap)  (4B) 

AD  =  ACd  (h  pV2  Ap)  (5B) 


5)  Calculate  the  corrected  lift,  CORL,  and  corrected  drag,  CORD,  response. 


CORL  =  AMP1  (channel  7)  -  AL  (6B) 

CORD  =  AMP1  (channel  12)  -  AD  (7E) 

The  above  procedure  was  repeated  for  the  second  harmonic  and  printed  in 
a  table  on  the  last  3  lines  of  page  2  output  (see  Fig.  2B) .  The  only 
significant  correction  to  the  data  occurs  to  the  lift  response  at  the 
fundamental  frequency  and  even  that  is  very  small,  normally  2-3%  but  not 
more  than  6%.  All  other  corrections  are  too  small,  in  the  third  significant 
digit,  to  be  meaningful.  Note  that  no  correction  is  applied  to  the  phase  angle 
of  the  response. 

The  correction  procedure  outlined  above  was  applied  in  two  ways.  First 
an  average  lift  amplitude  correction  was  applied  to  the  uncorrected  data 
representing  the  entire  run.  This  lift  correction  appears  in  channel  9  on 
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page  2  and  the  corrected  lift  amplitudes  are  in  channel  10.  Next,  a  point  by  point 
correction  was  performed.  Twenty  data  spots  were  taken  each  second,  and  each  data 
spot  was  individually  corrected.  Both  the  lift  and  drag  amplitudes  were  corrected 
in  this  manner  and  these  corrections  appear  in  the  table  on  the  last  three  lines  of 
page  2.  Throughout  the  test,  only  the  fundamental  frequency  lift  amplitude 
correction  was  at  all  significant,  all  other  corrections  being  too  small.  The 
difference  between  the  spot  by  spot  correction  method  and  average  correction  for 
the  run  method  is  infinitesimal  and  is  due  to  computer  round  off  error.  For  the 
sample  case  it  was  83.10  pounds  (369.6  N)  AMP1  lift  for  the  average  run  method  vs 
83.12  pounds  (369.7  N)  for  the  spot  by  spot  method. 


A  line  by  line  description  of  the  output  follows: 

Page  1  ~  Fig*  ^ 

Line  1  Experiment  title  &  date. 

Line  2  &  3  Title  identifying  the  run. 

Line  4-7  Shows  the  foil  depth-chord  ratio  measured  to  the  \  chord  location  of 
the  foil,  the  kinematic  viscosity,  mass  density  of  water  and  the 
water  depth. 


Line  8 


Lists  encounter  frequency  FE  in  hertz,  circular  encounter  frequency 
WE  in  radians/sec.  and  the  encounter  period  TE  in  seconds. 


Line  9 


Line  10 


Line  11 


Line  12 


Lists  the  wave  frequency  FO  in  hertz,  wave  circular  frequency  WO  in 
rad. /sec.  and  wave  period  TO  in  seconds. 

Lists  the  wave  length  LAMDA  in  feet  (0.3048  m) ,  the  wave  slope  in 
degrees  and  the  wave  celerity  in  knots. 

Lists  the  wave  length  to  foil  chord  ratio  LAMDA/L,  the  wave  steepness 
and  the  wave  length  to  wave  height  ratio,  LAMDA/H 

Lists  the  heading  (180°  indicates  head  seas) 


Line  13  &  14  Lists  the  actual  measured  speed  in  knots  and  the  nominal  speed  that 
was  specified  for  the  run. 


Line  15 


Lists  the  run  number  and  the  duration  of  the  data  collection  time  in 
seconds. 


Line  16  -  28  Lists  a  table  the  heading  of  which  is 

CHAN  -  data  channel 

CALIB  -  calibration  factor 

GAIN  -  gain  factor 

MEAN  -  algebraic  mean  of  the  data 
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The  individual 
Channel  1 
Channel  2 
Channel  3 
Channel  4 

Channel  5 
Channel  6 
Channel  7 
Channel  8 
Channel  9 
Channel  10 
Channel  11 
Channel  12 

Line  29 

Line  30 

Line  31 


Line  1-  15 
Line  16  -  18 

Line  19  -  29 


STDEN  -  standard  deviation  of  the  data 
R00TQ0  -  standard  deviation  multiplied  by  /2 . 0 

channels  of  interests  are: 

Drag  measured  in  pounds  (4.448  N) 

Lift  aft  measured  in  pounds  (4.448  N) 

Lift  fwd  measured  in  pounds  (4.448  N) 

Wave  height  measured  in  inches  (2.54  cm)  Wave  height  was  taken 
13.54  ft.  (4.126  m)  ahead  of  the  mid-chord  hydrofoil  position  on 
the  centerline 

Spare  wave  height  channel  -  not  used 
Carriage  speed  -  knots 

Uncorrected  lift  equal  to  the  sum  of  channels  2  and  3 
Not  used 

Lift  correction  due  to  sting  deflection 
Corrected  lift  equal  to  channel  7  minus  channel  9 
Not  used 

Drag  corrected  for  dynamometer  interaction 


Lists  KE  -  the  reduced  frequency  of  encounter,  FNC  -  the  chord  Froude 
Number,  FNH  -  the  foil  depth  Froude  number,  RN  -  the  Reynold's 
number  based  on  foil  chord. 

Lists  CF  -  the  Schoemerr  friction  coefficient,  CL  (script  symbol 
C^)  -  the  lift  coefficient  based  on  the  mean  value  of  channel  10, 
and  CD  -  (script  symbol  C  )  the  drag  coefficient  based  on  the  mean 
value  of  channel  12 

Lists  VFPS,  -  the  carriage  speed  in  feet/seconds  (0.3048  m/s)  and 
DAS  in  degrees.  DAS  'is  the  calculated  incremental  mean  angle  of 
attack  due  to  the  angular  deflection  of  the  dynamometer  nose  piece 
resulting  from  the  applied  mean  lift  as  recorded  in  channel  10. 

Page  2  -  Figure  2B  -  Unsteady  output  for  runs  in  waves  only 
Repeat  the  output  on  page  1 

Lists  the  harmonic  frequency  multiple  that  was  analyzed,  the 
harmonic  frequencies  in  rad. /sec.  and  the  number  of  cycles  of  data 
which  were  analyzed. 

Consists  of  a  table.  The  table  headings  are: 

R00TQ0  -  equal  to  the  standard  deviation  multiplied  by  /2 . 0 
AMP1/RQ0  -  is  the  quotient  of  the  amplitude  of  the  fundamental 
frequency  response  and  R00TQ0 

TRN  FUN  -  Truncation  Function,  disregard  this  quantity 
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AMP1  &  PHASE  -  amplitude  and  phase  at  the  fundamental  frequency  as 
determined  by  harmonic  analysis  for  each  channel  of 
data 

AMP2  &  PHASE  -  amplitude  and  phase  at  twice  the  fundamental  frequency 
as  determined  by  harmonic  analysis  for  each  data 
channel . 

The  data  channels  are  the  same  as  on  page  1  lines  17  -  28. 

Line  30  -  32  Starting  with  KE  these  are  the  same  as  on  page  1  lines  29  -  31. 

Line  33  CLR  -  fundamental  frequency  lift  response  operator  script  symbol; 

^L1  2 

-  =  (L /^pA  V  )  '  (1.0/ r)  where  L.  is  the  channel  10  AMPl 

rip  1 

value  as  it  appears  on  line  27. 


CDR  -  fundamental  frequency  drag  response  operator  script  symbol; 

*^D1  2 

-  =  (D1  /%  V  A  )  ‘  (1.0/r)  where  D  is  the  channel  12  AMPl 

rip  1 

value  as  it  appears  on  line  29 

Line  34  -  36  Consists  of  a  snort  table  of  additional  reduced  data  at  the 

fundamental  and  second  harmonic  frequency.  Quantities  in  the  lift 
and  drag  response  correction  due  to  foil  incidence  angle  deflections 
are  printed.  All  the  quantities  in  this  table  were  obtained  on  a 
point  basis  as  opposed  to  the  average  lift  and  drag  corrections 
shown  on  channels  9,  10  &  12. 

The  table  headings  are: 

HARM  -  harmonic  of  the  fundamental  frequency 

DA  -  the  calculated  amplitude  of  the  angular  deflection  of  the  nose 
piece  (in  degrees)  due  to  the  oscillating  lif'  forces  imposed 
on  the  foil.  It  is  calculated  by  taking  static  angular 
deflection  calibration  of  the  nose  piece  under  lift  and 
multiplying  by  the  oscillating  lift  force  (AMPl  channel  10  for 
the  fundamental  frequency  or  AMP2  for  the  second  harmonic) 

DCL  -  the  change  in  lift  coefficient  calculated  by  multiplying  DA 
with  the  foil  lift  curve  slope  obtained  from  previously 
unpublished  smooth  water  experiments  on  this  foil 

2 

DCD  -  the  change  in  drag  coefficient  calculated  by  multiplying  (DA) 
by  the  3CD  obtained  from  previously  unpublihsed  smooth  water 


experiments . 

DL  -  oscillating  lift  amplitude  correction  based  on  DCL 
DD  -  oscillating  drag  amplitude  correction  based  on  DCD 
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CORL  -  Corrected  value  of  the  oscillating  lift  amplitude  obtained 
by  subtracting  DL  from  the  channel  7  amplitude  -  the  CORL 
value  should  be  almost  identical  to  the  channel  10  value  of 
the  corrected  lift. 

CORD  -  Corrected  value  of  the  oscillating  drag  amplitude 

CL  -  Amplitude  of  the  oscillating  lift  in  coefficient  form  based  on 
CORL 

D  -  Amplitude  of  the  oscillating  drag  in  coefficient  form  based  onCORD 

CLR  -  Lift  response  operator  -  lift  based  on  CL  line  34 

CDR  -  Drag  response  operator  -  drag  based  on  CD  line  34 


i 
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OTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INIEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


